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We predict analytically and confirm with numerical simulations that inter-mode dispersion in 
nanowire waveguide arrays can be tailored through periodic waveguide bending, facilitating flexible 
spatio-temporal reshaping without break-up of femtosecond pulses. This approach allows simulta¬ 
neous and independent control of temporal dispersion and spatial diffraction that are often strongly 
connected in nanophotonic structures. 


High-index-contrast nanowires offer unique advantages 
for manipulation of optical pulses in compact photonic 
circuits, providing high field confinement and enabling 
precise dispersion engineering mm- In particular op¬ 
tical chips based on silicon sub-wavelength waveguides 
allow for efficient frequency conversion mm, all-optical 
pulse control [5] and all-optical switching [6]. Further¬ 
more, couplers [ 7 ] and arrays of coupled nanowire waveg¬ 
uides [SIE] open possibilities for efficient spatio-temporal 
shaping of optical pulses. In order to harness these oppor¬ 
tunities, it is essential to develop approaches to simulta¬ 
neously and independently control temporal and spatial 
dispersion, as these characteristics can be strongly con¬ 
nected in nanophotonic structures. This connectivity can 
lead to difficulties designing a waveguide array support¬ 
ing propagation of ultrashort pulses, since pulses either 
disperse due to strong temporal dispersion or break-up 
due to strong spatial diffraction [9]. 

One possible approach to achieve required spatio- 
temporal dispersion is to carefully design waveguide ar¬ 
ray geometry and use complex photonic crystal struc¬ 
tures m- That is however a very complicated method. 
Another way to tailor dispersion is by introducing peri¬ 
odic waveguide bending [IIHIS]. This approach allows 
relatively simple fabrication and offers substantial design 
flexibility. Periodic waveguide bending was introduced 
as an effective tool for polychromatic diffraction man¬ 
agement UBIIHlIii], however it has only been studied 
in the context of continuous light illumination and con¬ 
ventional micro-scale waveguides. In this work, we de¬ 
velop an approach to simultaneously control spatial and 
temporal dispersion and demonstrate through numerical 
simulations the application of this concept to the sup¬ 
pression of ultrashort pulse distortion and break-up in 
nano-waveguide arrays. 

Ultrashort pulses have a broad spectrum encompassing 
a large range of frequencies cj, therefore temporal disper¬ 
sion characterized by a propagation constant l3s{oo) has a 
significant influence on the pulse dynamics. In waveguide 
arrays pulses can also switch between different waveg¬ 
uides. One of the most important parameters charac¬ 
terizing waveguide arrays is a coupling coefficient Cs{uj), 
which determines the rate at which light couples between 
the neighboring waveguides and thus regulates the spatial 


dispersion. The coupling coefficient for straight lossless 
waveguides is real: Cs{uj) = Re[Cs(co’)]. 

It was shown previously m that bending the waveg¬ 
uides affects the coupling by introducing an additional 
phase shift. If all waveguides in an array have the same 
bending profile ^ 0 ( 2 :), where xq is a transverse coordinate 
of the waveguide center and z is the propagation direc¬ 
tion, then the complex electric field amplitude in n-th 
waveguide of the array satisfies the following coupled¬ 
mode equations [13]: 

ldzEn{z,0j) + l3s{uj)En{z,0j) = 

- C's(w)exp [inQdwXQ{z)u!/c]En-i{z,u!) (1) 

- C's(w)exp [-tnodwXo{z)u/c]En+i{z,Lo). 

Here no is an effective refractive index, ui = ^-KcjX is the 
angular frequency, A is a light wavelength in vacuum, and 
dw is the distance between the coupled waveguides. 

We extend this method to consider the dynamics of 
ultrashort pulses and study its applicability to nanopho¬ 
tonic structures. In low-index waveguide arrays both 
the propagation constant of individual waveguides and 
the coefficient characterizing coupling between the neigh¬ 
boring waveguides are mildly dispersive. In contrast, 
in nanowire high-index waveguides the dispersion can 
be much stronger, and also small changes of waveguide 
cross-section dramatically affect both temporal disper¬ 
sion and spatial diffraction [9|. 

To investigate the pulse dynamics in the coupled 
nanowires we combine the approaches previously devel¬ 
oped for the description of nanowire arrays [9] and curved 
conventional waveguide arrays [Hill]. We derive the fol¬ 
lowing system of equations by applying to Eq. 0 the 
Fourier transform En{z,t) = J dujEn{z,uj) ex.p {—luot), 
where t is time, and perform a Taylor expansion of cou¬ 
pling and propagation coefficients: 

id^En{z,t) + pEn{z,t) = 

- C{z)En-l{z,t) - C\z)En+l{z,t). 

Here $ determines the temporal dispersion in a waveg¬ 
uide, and C characterizes the coupling between the neigh- 
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FIG. 1. (a) The scheme of three straight coupled nanowires, 

(b) Coupling coefficient vs. wavelength between the neigh¬ 
boring waveguides, (c) Supermode propagation constants vs. 
wavelength (solid lines), and the corresponding profiles of the 
supermodes (cross sections of three coupled waveguides). 
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FIG. 2. Pulse intensity evolution along straight coupled 
nanowires: (a) left, (b) central, and (c) right nanowire. The 
temporal axis corresponds to a moving time frame with the 
group velocity at the central wavelength, r — t — z/Si . 


boring waveguides: 

M Q M / \ 

/3 = E C{z) = ^ (3) 

ml ml 

m=0 m=0 

where M is a sufficiently large number to capture the 
dispersion features over the pulse bandwidth. The Taylor 
coefficients are 

Cm{z) = [Cs{uj) , 

where ujq is the central pulse frequency. Note that the dis¬ 
persion coefficients {Prn) do not depend on the coupling. 
On the other hand, the coupling coefficients (c^) depend 
non-trivially on the propagation distance {z) through an 
interplay between the dispersion of coupling between the 
straight waveguides (C^) and the bending profile [^ 0 ( 2 :)] 
induced dispersion. 

It was shown m that for polychromatic light propa¬ 
gation in periodically curved waveguide arrays, after the 
full bending period the beam diffraction is the same as 
in a straight array with the effective coupling coefficient. 
We check that the same approach can be applied to the 
pulse propagation when ^ 0 ( 2 :) = xo{z L^) and I /5 is a 
modulation period, and the effective coupling is: 

r^b 

= Cs{uj)L^^ / cos[nodwXo{z)uj/c]dz. (5) 
Jo 

The corresponding Taylor expansion coefficients are: 

cf = ( 6 ) 

We see that diffraction of beams is defined by an inter¬ 
play of the additional bending-induced dispersion intro¬ 
duced through the frequency dependence of the integral 
in Eq. and the intrinsic frequency dependence of the 
coupling coefficient in a straight waveguide array Cs{uj). 


We investigate the influence of the periodic waveguide 
bending on the pulse reshaping, and consider a represen¬ 
tative example of a cosine profile with the amplitude A 
and period Li^: 

xo{z) = Acos{27rz/Li)). (7) 

Below we show that a special combination of the mod¬ 
ulation parameters allows us to suppress the dispersion 
of the effective coupling coefficient and accordingly avoid 
the pulse distortion. 

To demonstrate our approach, we consider coupled Si 
nanowire waveguides and use COMSOL RF module to 
calculate the dispersion in straight coupled waveguides. 
The dimensions of the wires are as follows. The wires 
are 220 nm high and 330 nm wide, placed on a silica slab. 
There is a 100 nm high etching mask with refractive index 
1.35 on top of wires. Otherwise the wires are surrounded 
by air. We choose these parameters to obtain nearly zero 
group velocity dispersion ^^2 ~ 0 in the proximity of Aq = 
1.5/im wavelength for a single nanowire, as this would 
minimize the pulse distortion. We however emphasize 
that even in this regime, a pulse can exhibit distortion 
due to the dispersion of coupling between the nanowires 
in an array [9]. 

To determine the coupling strength between the neigh¬ 
boring waveguides, we follow the approach of Ref. [9] and 
analyze a two-waveguide coupler with wire-to-wire sep¬ 
aration of 330 nm. We calculate the propagation con¬ 
stants for symmetric and antisymmetric supermodes of 
the coupler and respectively. The 

propagation constant for a single waveguide can be well 
approximated by the average of the symmetric and an¬ 
tisymmetric supermode propagation constants, /3s{uj) ^ 
j^sym(^) while the difference defines the 

coupling coefficient C,(cj) = The 

coefficients of the Taylor expansion of the propagation 
constant are found as /3o = 7.85/im“^, = 17fs/im“^, 

(32 = 1.4fs^/im“^. Accordingly, we find that the group 
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FIG. 3. (a) The scheme of three periodically bent coupled 
nanowires with the bending period Lb = Lc = 15.75 
(b) Effective coupling coefficient over one bending period be¬ 
tween the neighboring nanowires vs. wavelength, (c) Average 
supermode propagation constants over one bending period vs. 
wavelength. 
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FIG. 4. Pulse intensity evolution along periodically curved 
coupled nanowires: left edge nanowire (a), central nanowire 
(b), right edge nanowire (c). The temporal axis corresponds 
to a moving time frame with the group velocity at the central 
wavelength, r — t — z/3i . 


velocity dispersion is indeed suppressed for pulses with 
duration down to 100 fs. 

Next we consider an array of three coupled nanowires 
with all parameters as noted before, see Fig. 0a). We 
use a three-waveguide system for all the following calcu¬ 
lations, since it allows us to fully demonstrate the pulse 
control method proposed in this work. In agreement 
with predictions of Ref. [9], we notice that the varia¬ 
tions of the coupling coefficient for such waveguide ar¬ 
ray are significant across a relatively narrow spectrum 
[see Fig. Bb)], which could lead to temporal reshaping of 
short pulses during propagation. The Taylor expansion 
of the coupling coefficient for straight waveguides is cq = 
69.7 mm“^, ci = —0.71 fs /im“^, C 2 = 3.6 fs^ /im“^, which 
reveals strong linear dispersion and small quadratic dis¬ 
persion in the wavelength range between 1.46/rm and 
1.54/im. The coupled waveguides support three super¬ 
modes, which characteristic spatial profiles are shown in 
Fig. 0 c). We calculate the dependence of the supermode 
propagation constants on the wavelength, see Fig. [^c). 
These dependencies have different slopes corresponding 
to different supermode velocities. We show below that 
this leads to pulse splitting, which can be suppressed via 
periodic waveguide bending. 

First we analyze the pulse dynamics in straight waveg¬ 
uides with length L = 20Lc = 315 /im, where Lc = 
15.75 /im is the length required for full coupling from one 
waveguide to another at the wavelength Aq. As an input, 
we consider a single 100 fs long transform-limited Gaus¬ 
sian pulse with the central wavelength Aq coupled to the 
left nanowire of the straight waveguide array. Figs, [^a- 
c) demonstrate that initially the pulse couples from the 
left (a) to the central (b) and then to right (c) nanowire 
without significant distortions. Then the pulse starts to 
split into three separate pulses in the edge waveguides 
(a,c) and into two pulses in the central waveguide (b), 
in agreement with the previous study [9]. These pulses 


propagate with different group velocities, which corre¬ 
spond to three different supermode velocities supported 
in the structure [see Fig. [^c)]. Such behavior demon¬ 
strates that although the single nanowire propagation 
dispersion can be engineered, spatial diffraction in ar¬ 
rays of nanowires is still strongly limited by the coupling 
dispersion. Moreover the propagation dispersion and the 
coupling dispersion in nanowire waveguide arrays are in¬ 
terconnected, and therefore an approach allowing for the 
independent control of these characteristics would offer 
essential benefits for various applications. 

Next, we investigate the influence of the periodic 
waveguide bending on the pulse reshaping. We choose 
a bending profile according to Eq. 0 . We vary the 
bending amplitude A and search for the minima of the 
coupling dispersion dCeff/duj in the vicinity of Aq. We 
choose the bending period L^ = Lc = 15.75 /im, as it al¬ 
lows us to consider nanowires with smaller curvature for 
the purposes of easier potential fabrication and reduc¬ 
tion of propagation losses. As we show bellow, one can 
choose a bending profile that simultaneously allows for a 
strong coupling dispersion control and does not introduce 
bending propagation losses. 

We calculate the effective coupling coefficient Cef / us¬ 
ing Eq. 0- We choose the value of A = Amin corre¬ 
sponding to the first minimum of dCef f /du, and accord¬ 
ingly the smallest suitable bending curvature. The opti¬ 
mal bending amplitude is found to be equal to Amin = 
1.3/im. The bending losses for the corresponding cur¬ 
vature value should be practically absent according to 
the previous studies of bent nanowire waveguides m- 
The resulting effective coupling coefficient Ceff shown 
in Fig. I^b) becomes almost constant over a broad spec¬ 
tral region in comparison to that for the straight waveg¬ 
uides [see Fig. 0 b)]. The Taylor expansion of the ef¬ 
fective coupling coefficient is = 9.96 mm“^, cf^ = 
—0.0076 fs/im“^, = —14.2 fs^/im“^. Although the 
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quadratic coupling dispersion is slightly increased in com¬ 
parison to the straight waveguide array, the linear cou¬ 
pling dispersion, which is the main temporal reshaping 
driver for 100 fs long pulses in such structures, is sup¬ 
pressed by two orders of magnitude compared to the 
straight waveguides. In Fig. |^c) we plot the supermode 
propagation constants for the curved waveguide arrays 
with three coupled nanowires calculated with the use of 
the effective coupling coefficient. The propagation con¬ 
stants for the three supermodes now have similar slopes, 
which suggests that short-pulse break-up due to coupling 
would be suppressed. 

We now calculate the intensity evolution of a 100 fs 
transform-limited pulse coupled to the left nanowire 
of the periodically curved waveguide array using the 
Eqs. ([^-Q. We see in Figs, [^a-c) that as a result of 
vanishing coupling dispersion the temporal pulse break¬ 
up is suppressed, and a pulse can now be switched as a 
whole between the waveguides. Thus the temporal and 
the spatial dispersion in nanowire waveguide arrays can 
be controlled independently via single waveguide disper¬ 
sion engineering and periodic waveguide bending. 

These results demonstrate that spatio-temporal dis¬ 
persion engineering in high-index-contrast nanowire 
waveguide arrays can be efficiently realized through the 
introduction of periodic waveguide bending, which can 
enable flexible spatio-temporal manipulation of femtosec¬ 
ond pulses. We anticipate that these results will open 
novel approaches to on-chip all-optical light control [8]. 
This approach can also be useful for enhanced paramet¬ 
ric frequency conversion [2] and broadband photon-pair 
generation and quantum walks [16]. 
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